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THE MILLISECOND MAGNETAR CENTRAL ENGINE IN SHORT GRBS 
Hou-Jun Lu\ Bing Zhang\ Wei-Hua Lei^, Ye LL, Paul D Lasky^’'* 

ABSTRACT 

One favored progenitor model for short duration gamma-ray bursts (GRBs) is the coalescence of 
two neutron stars (NS—NS). One possible outcome of such a merger would be a rapidly spinning, 
strongly magnetized neutron star (known as a millisecond magnetar). These magnetars may be 
“supra-massive,” implying that would collapse to black holes after losing centrifugal support due to 
magnetic dipole spin down. By systematically analyzing the Burst Alert Telescope (BAT)-XRT light 
curves of all short GRBs detected by Swift, we test how consistent the data are with this central engine 
model of short GRBs. We find that the so-called “extended emission” feature observed with BAT in 
some short GRBs is fundamentally the same component as the “internal X-ray plateau” as observed 
in many short GRBs, which is defined as a plateau in the light curve followed by a very rapid decay. 
Based on how likely a short GRB is to host a magnetar, we characterize the entire Swift short GRB 
sample into three categories: the “internal plateau” sample, the “external plateau” sample, and the 
“no plateau” sample. Based on the dipole spin down model, we derive the physical parameters of the 
putative magnetars and check whether these parameters are consistent with expectations from the 
magnetar central engine model. The derived magnetar surface magnetic field Bp and the initial spin 
period Pq fall into a reasonable range. No GRBs in the internal plateau sample have a total energy 
exceeding the maximum energy budget of a millisecond magnetar. Assuming that the beginning of 
the rapid fall phase at the end of the internal plateau is the collapse time of a supra-massive magnetar 
to a black hole, and applying the measured mass distribution of NS—NS systems in our Galaxy, we 
constrain the neutron star equation of state (EOS). The data suggest that the NS EOS is close to the 
GMl model, which has a maximum non-rotating NS mass of Mtov ~ 2.37Mq. 

Subject headings: gamma rays: general- methods: statistical- radiation mechanisms: non-thermal 


1. INTRODUCTION 

Gamma-ray bursts (GRBs) are classified into cate¬ 
gories of “long soft” and “short hard” (SGRB) based on 
the observed duration (Tgo) and hardness ratio of their 
prompt gamma-ray emission (Kouveliotou et al. 1993). 
Long GRBs are found to be associated with core-collapse 
supernovae (SNe; e.g. Galama et al. 1998; Hjorth et al. 
2003; Stanek et al. 2003; Campana et al. 2006; Xu et 
al. 2013), and typically occur in irregular galaxies with 
intense star formation (Fruchter et al. 2006). They are 
likely related to the deaths of massive stars, and the “col- 
lapsar” model has been widely accepted as the standard 
paradigm for long GRBs (Woosley 1993; MacFadyen & 
Woosley 1999). The leading central engine model is a 
hyper-accreting black hole (e.g. Popham et al. 1999; Lei 
et al. 2013). Alternatively, a rapidly spinning, strongly 
magnetized neutron star (millisecond magnetar) may be 
formed during core collapse. In this scenario, magnetic 
fields extract the rotation energy of the magnetar to 
power the GRB outflow (Usov 1992; Thompson 1994; 
Dai & Lu 1998; Wheeler et al. 2000; Zhang & Meszaros 
2001; Metzger et al. 2008, 2011; Lyons et al. 2010; Buc- 
dantini et al. 2012; Lii & Zhang 2014). 

In contrast, short GRBs are found to be associated 
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with nearby early-type galaxies with little star formation 
(Barthelmy et al. 2005; Berger et al. 2005; Gehrels et al. 
2005; Bloom et al. 2006), to have a large offset from the 
center of the host galaxy (e.g. Fox et al. 2005; Fong et 
al. 2010), and to have no evidence of an associated SNe 
(Kann et al. 2011, Berger 2014 and references therein). 
The evidence points toward an origin that does not in¬ 
volve a massive star. The leading scenarios include the 
merger of two neutron stars (NS—NS, Paczynski 1986; 
Eichler et al. 1989) or the merger of a neutron star and 
a black hole (Paczynski 1991). For NS—NS mergers, the 
traditional view is that a BH is formed promptly or fol¬ 
lowing a short delay of up to hundreds of milliseconds 
(e.g. Rosswog et al. 2003; Rezzolla et al. 2011; Liu et 
al. 2012). Observations of short GRBs with Swift, on 
the other hand, indicated the existence of extended cen¬ 
tral engine activity following at least some short GRBs 
in the form of extended emission (EE; Norris & Bonnel 
2006), X-ray flares (Barthelmy et al. 2005; Campana 
et al. 2006), and, more importantly, “internal plateaus” 
with rapid decay at the end of the plateaus (Rowlinson 
et al. 2010, 2013). These observations are difficult to 
interpret within the framework of a black hole central 
engine, but are consistent with a rapidly spinning mil¬ 
lisecond magnetar as the central engine (e.g. Dai et al. 
2006; Gao & Fan 2006; Metzger et al. 2008; Rowlinson 
et al. 2010, 2013; Gompertz et al. 2013, 2014). 

About 20% of short GRBs detected with Swift have 
EE (Sakamoto et al. 2011) following their initial short, 
hard spike. Such EE typically has a lower flux than the 
initial spike, but can last for tens of seconds (e.g. Perley 
et al. 2009). The first short GRB with EE detected with 
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Swift was GRB 050724, which had a hard spike of Tqo ~ 
3s followed by a soft tail with a duration of ~150 s in 
the Swift Burst Alert Telescope (BAT; Barthelmy et al. 
2005) band. The afterglow of this GRB lies at the out- 
skirt of an early-type galaxy at a redshift of 2=0.258. It 
is therefore a “smoking-gun” burst of the compact star 
merger population (Barthelmy et al. 2005; Berger et al. 
2005). GRB 060614 is a special case with a light curve 
characterized by a short/hard spike (with a duration ^ 
5s) followed by a series of soft gamma-ray pulses lasting 
^100 s. Observationally, it belongs to a long GRB with¬ 
out an associated SNe (with very deep upper limits of 
the SN light, e.g. Della Valle et al. 2006; Fynbo et al. 
2006; Gal-Yam et al. 2006). Some of its prompt emis¬ 
sion properties, on the other hand, are very similar to a 
short GRB (e.g. Gehrels et al. 2006). Through simu¬ 
lations, Zhang et al. (2007b) showed that if this burst 
were a factor of 8 less luminous, it would resemble GRB 
050724 and appear as a short GRB with EE. Norris & 
Bonnell (2006) found a small fraction of short GRBs in 
the BATSE catalog qualitatively similar to GRB 060614. 
It is interesting to ask the following two questions. Are 
short GRBs with EE different from those without EE? 
What is the physical origin of the EE? 

Swift observations of the X-ray afterglow of short 
GRBs, on the other hand, provide some interesting clues. 
A good fraction of Swift short GRBs exhibit an X-ray 
plateau followed by a very sharp drop with a temporal 
decay slope of more than 3. The first case was GRB 
090515 (Rowlinson et al. 2010). It showed a nearly flat 
plateau extending to over 180 s before rapidly falling off 
with a decay slope of a ^ 130 Such rapid decay cannot 
be accommodated by any external shock model, and so 
the entire X-ray plateau emission has to be attributed to 
the internal dissipation of a central engine wind. Such 
an “internal plateau” has previously been observed in 
some long GRBs (e.g. Troja et al. 2007; Lyons et al. 
2010), but they are also commonly observed in short 
GRBs (Rowlinson et al. 2013). These plateaus can be 
interpreted as they internal emission of a spinning-down 
magnetar which collapses into a black hole at the end of 
the plateau (Troja et al. 2007; Rowlinson et al. 2010; 
Zhang 2014). 

If magnetars are indeed operating in some short GRBs, 
then several questions emerge. What fraction of short 
GRBs have a millisecond magnetar central engine? What 
are the differences between short GRBs with EE and 
those that have an internal plateau but no EE? Is the 
total energy of the magnetar candidates consistent with 
the maximum rotation energy of the magnetars according 
to the theory? What are the physical parameters of the 
magnetar candidates derived from observational data? 
How can one use the data to constrain the equation of 
state (EOS) of neutron stars? 

This paper aims to address these interesting questions 
through a systematic analysis of both Swift/BAT and X- 
ray Telescope (XRT) data. The data reduction details 
and the criteria for sample selection are presented in §2. 
In §3, the observational properties of short GRBs and 
their afterglows are presented. In §4, the physical pa¬ 
rameters of the putative magnetars are derived and their 

® The convention iT oc is adopted throughout the pa¬ 
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statistical properties are presented. The implications on 
the NS EOS are discussed. The conclusions are drawn in 
§5 with some discussion. Throughout the paper, a con¬ 
cordance cosmology with parameters = 71 km s“^ 
Mpc VLm = 0.30, and Ha = 0.70 is adopted. 

2. DATA REDUCTION AND SAMPLE SELECTION 
CRITERIA 

The Swift BAT and XRT data are downloaded from 
the Swift data archive0 We systematically process the 
BAT and XRT GRB data to extract light curves and 
time-resolved spectra. We developed an IDL script to au¬ 
tomatically download and maintain all of the Swift BAT 
data. The HEAsoft package version 6.1 0, including bate- 
convert, batbinevt Xspec, Xselect, Ximage, and the Swift 
data analysis tools are used for data reduction. The de¬ 
tails of the data analysis method can be found in several 
previous papers (Liang et al. 2007; Zhang et al. 2007c; 
Lii et al. 2014) in our group, and Sakamoto et al. (2008). 

We analyze 84 short GRBs observed with Swift be¬ 
tween 2005 January and 2014 August. Among these, 44 
short GRBs are either too faint to be detected in the 
X-ray band, or do not have enough photons to extract a 
reasonable X-ray light curve. Our sample therefore only 
comprises 40 short GRBs, including 8 with EE. 

We extrapolate the BAT (15-150 keV) data to the XRT 
band (0.3-10 KeV) by assuming a single power-law spec¬ 
trum (see also O’Brien et al. 2006; Willingale et al. 2007; 
Evans et al. 2009). We then perform a temporal fit to 
the light curve with a smooth broken power law in the 
rest frame 0 



to identify a possible plateau in the light curve. Here, 
tb is the break time. Ft = Fq ■ is the flux at the 

break time tb, ai and a 2 are decay indices before and af¬ 
ter the break, respectively, and w describes the sharpness 
of the break. The larger the u parameter, the sharper 
the break. An IDL routine named “mpfitfun.pro” is em¬ 
ployed for our fitting (Markwardt 2009). This routine 
performs a Levenberg-Marquardt least-square fit to the 
data for a given model to optimize the model parameters. 

Since the magnetar signature typically invokes a 
plateau phase followed by a steeper decay (Zhang & 
Meszaros 2001), we search for such a signature to de¬ 
cide how likely a GRB is to be powered by a magnetar. 
Similar to our earlier work (Lii & Zhang 2014), we intro¬ 
duce three grades to define the likelihood of a magnetar 
engine. 

• The internal plateau(Internal) sample: this sample 
is defined by those bursts that exhibit a plateau 

® http://www.swift.ac.uk/archive/obs.php?burst=l 

^ Another empirical model to fit GRB X-ray afterglow light 
curves is that introduced by Willingale et al. (2007, 2010). The 
function was found to be a good fit of the external plateaus of 
long GRBs (e.g. Dainotti et al. 2010), but cannot fit the internal 
plateaus which are likely due to the magnetar origin (e.g. Lyons et 
al. 2010). We have tried to use the Willingale function to fit the 
data in our sample, but the fits are not good. This is because our 
short GRB sample includes a large fraction of internal plateaus. 
We therefore do not use the Willingale function to fit the light 
curves in this paper. 
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followed by a decay with or steeper than 3. 
The decay is expected by the magnetar dipole 
spin down model (Zhang & Meszaros 2001), while 
a slope steeper than three is an indication that the 
emission is powered by internal dissipation of the 
magnetar wind, since essentially no external shock 
model can account for such a steep decay. This 
sample is similar to the “Gold” sample defined by 
Lii & Zhang (2014)0, but with the inclusion of 
two GRBs with a decay following the plateau. 
These two GRBs (GRB 061201 and GRB 070714B) 
also have a plateau index close to 0 as demanded by 
the magnetar spin down model, and therefore are 
strong candidates for magnetar internal emission. 
For those cases with a post-plateau decay index 
steeper than three, the rapid decay at the end of 
plateau may mark the implosion of the magnetar 
into a black hole (Troja et al. 2007; Zhang 2014). 
Altogether there are 20 short GRBs identified as 
having such a behavior, 13 of which have redshift 
measurements and 7 of which are short GRBs with 
EE. For these latter GRBs, the extrapolated X-ray 
light curves from the BAT band in the EE phase 
resemble the internal plateaus directly detected in 
the XRT band in other GRBs. The light curves of 
these 22 GRBs are presented in Fig.l, along with 
the smooth broken power-law fits. The fitting pa¬ 
rameters are summarized in Table 1. 

• The external plateau (External) sample: this sam¬ 
ple includes those GRBs with a plateau phase fol¬ 
lowed by a normal decay segment with a post-decay 
index close to -1. The pre- and post-break temporal 
and spectral properties are consistent with the ex¬ 
ternal forward shock model with the plateau phase 
due to continuous energy injection into the blast- 
wave. This sample is similar to the Silver and Alu¬ 
minum samples in Lii & Zhang (2014). We identi¬ 
fied 10 GRBs in this group 0 The XRT light curves 
are presented in Figure 2 along with the smooth 
broken power-law fits. The fitting results are pre¬ 
sented in Table 1. 

• No plateau (Non) sample: we identify 8 GRBs that 
do not have a significant plateau behavior. They 
either have a single power-law decay, or have er¬ 
ratic flares that do not present a clear magnetar 
signature. 

Figure 3 collects all of the light curves of the GRBs in 
our samples. The Internal sample with or without EE is 
collected in Figures 3(a) and 3(b); the External sample 
(without EE) is collected in Figure 3(c); and the Non 
sample are collected in Figure 3(d). 

® Lii & Zhang (2014) studied the magnetar engine candidates 
for long GRBs. The grades defined in that paper were based on 
the following criteria. Gold sample: those GRBs which dispaly an 
“internal plateau”; Silver sample: those GRBs which display an 
“external plateau,” whose energy injection parameter q is consis¬ 
tent with being 0, as predicted by the dipole spin down model of 
GRBs; Aluminum sample: those GRBs which display an external 
plateau, but the derived q parameter is not consistent with 0; Non- 
magnetar sample: those GRBs which do not show a clear plateau 
feature. 

® The SN-less long GRB 060614 is included in this category. It 
has EE and an additional external plateau at late times. 


3. DERIVED PHYSICAL PARAMETERS AND STATISTICS 

In this section, we derive physical parameters of the 
short GRBs in various samples, and perform some statis¬ 
tics to compare among different samples. 

3.1. EE and internal plateau 

Our first task is to investigate whether short GRBs 
with EE are fundamentally different from those without 
EE. The EE has been interpreted within the magnetar 
model as the epoch of tapping the spin energy of the 
magnetar (Metzger et al. 2008; Bucciantini et al. 2012). 
On the other hand, a good fraction of short GRBs with¬ 
out EE have an internal plateau that lasts for hundreds 
of seconds, which can also be interpreted as the inter¬ 
nal emission of a magnetar during the spin-down phase 
(Troja et al. 2007; Yu et al. 2010; Rowlinson et al. 
2013; Zhang 2014). It would be interesting to investi¬ 
gate whether or not there is a connection between the 
two groups of bursts. 

Analyzing the whole sample, we find that the short 
GRBs with EE do not show a plateau in the XRT band 
(except GRB 060614, which shows an external plateau 
at a later epoch). Extrapolating the BAT data to the 
XRT band, the EE appears as an internal plateau (Fig¬ 
ure. 1). Fitting the joint light curve with a broken power- 
law model, one finds that there is no significant differ¬ 
ence in the distribution and cumulative distribution of 
the plateau durations for the samples with and with¬ 
out EE (Figure 4(a)). The probability (p) that the two 
samples are consistent with one another, as calculated 
using a student’s t-test, is 0.6513 Figure 4(b) shows the 
redshift distributions of those short GRBs in our sam¬ 
ple which have redshift measurements. Separating the 
sample into EE and non-EE sub-samples does not re¬ 
veal a noticeable difference. In Figure 4(c), we show 
the flux distribution of the plateau at the break time. 
It is shown that the distribution for the EE sub-sample 
(mean flux logF), = —8.74± 0.12 ergs s“^ cm“^) is sys¬ 
tematically higher than that for the non-EE sub-sample 
(mean flux logF;, = —9.84 ± 0.07 ergs s“^ cm“^). How¬ 
ever, the combined sample (Figure 4(d)) shows a single¬ 
component log-normal distribution with a mean flux of 
logFh = —9.34 ± 0.07 ergs s“^ cm“^, with a student’s 
t-test probability of p = 0.76 of belonging to the same 
parent sample. This suggests that the EE GRBs are 
simply those with brighter plateaus, and the detection of 
EE is an instrumental selection effect. We also calculate 
the luminosity of the internal plateau at the break time 
for both the GRBs with and without EE. If no redshift 
is measured, then we adopt z =0.58, the center value 
for the measured redshift distribution (Figure 4(b)). We 
find that the plateau luminosity of the EE (logLo = 
49.41 ± 0.07 ergs s“^) is systematically higher than the 
no-EE sample (logLo = 48.68±0.04 ergs s“^), see Figure 
4(e). However, the joint sample is again consistent with a 
single component (logLo = 48.91 ±0.07 ergs s“^, Figure 
4(f)), with a student’s t-test probability of p = 0.74. For 
the samples with the measured redshifts only, our results 
(shown in the inset of Figure 4(e) and 4(f)), the results 
are similar. 

The hypothesis that the two distributions are from a same 
parent sample is statistically rejected if p < 0.05. The two samples 
are believed to have no significant difference if p > 0.05. 
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The distributions of the plateau duration, flux, and lu¬ 
minosity suggest that the EE and X-ray internal plateaus 
are intrinsically the same phenomenon. The different 
plateau luminosity distribution, along with the similar 
plateau duration distribution, suggest that the fraction 
of short GRBs with EE should increase with softer, more 
sensitive detectors. The so-called “EE” detected in the 
BAT band is simply the internal plateau emission when 
the emission is bright and hard enough. 

3 . 2 . The host offset and local environment of Internal 
and External samples 

One curious question is why most (22) short GRBs 
have an internal plateau, whereas some others (10) show 
an external plateau. One naive expectation is that the 
External sample may have a higher circumburst density 
than the Internal sample, and so the external shock emis¬ 
sion will be greatly enhanced. It has been found that 
short GRBs typically have a large offset from their host 
galaxies (Fong et al. 2010; Fong & Berger 2013; Berger 
2014), so that the local interstellar medium (ISM) den¬ 
sity may be much lower than that of long GRBs (e.g. 
Fan et al. 2005; Zhang et al. 2009; Kann et al. 2011). 
This is likely due to the asymmetric kicks during the 
SNe explosions of the binary systems when the two com¬ 
pact objects (NS or BH) were born (e.g. Bloom et al. 
1999, 2002). If the circumburst density is the key factor 
in creating a difference between the Internal and Exter¬ 
nal samples, then one would expect that the offset from 
the host galaxy is systematically smaller for the External 
sample than the Internal sample. 

With the data collected from the literature (Fong et 
al 2010, Leibler & Berger 2010, Fong & Berger 2013, 
Berger 2014), we examine the environmental effect of 
short GRBs within the Internal and External samples. 
The masses, ages, and specific star formation rates of the 
host galaxies do not show statistical differences between 
the two samples. The physical offsets and the normal¬ 
ized offset43 of these two samples are shown in the left 
and right panels of Figure 5. It appears that the ob¬ 
jects in the External sample tend to have smaller offsets 
than those in the Internal sample, both for the physical 
and normalized offsets. This is consistent with the above 
theoretical expectation. Nonetheless, the two samples 
are not well separated in the offset distributions. Some 
GRBs in the External sample still have a large offset. 
This may suggest a large local density in the ISM or 
intergalactic medium (IGM) far away from the galactic 
center, or that some internal emission of the nascent mag- 
netars may have observational signatures similar to the 
external shock emission. 

3 . 3 . Energetics and luminosity 

Similar to Lii & Zhang (2014), we derive the isotropic 
7 -ray energy {E^^iso) and isotropic afterglow kinetic en¬ 
ergy (EK,iso) of all of the short GRBs in our sample. To 
calculate E^^iso, we use the observed fluence in the de¬ 
tector’s energy band and extrapolate it to the rest-frame 
1 — 10"^ keV using spectral parameters with fc-correction 

The normalized offsets are defined as the physical offsets nor¬ 
malized to re, the characteristic size of a galaxy defined by Equa¬ 
tion (1) of Fong et al. (2010). 


(for details, see Lii & Zhang, 2014). If no redshift is 
measured, then we use z =0.58 (see Table 2). 

To calculate Ekjso, we apply the method described in 
Zhang et al. (2007a). Since no stellar wind environment 
is expected for short GRBs, we apply a constant density 
model. One important step is to identify the external 
shock component. If an external plateau is identified, 
then it is straightforward to use the afterglow flux to de¬ 
rive Ek^iso- The derived Ekjso is constant during the 
normal decay phase, but it depends on the time during 
the shallow decay phase (Zhang et al. 2007a). We there¬ 
fore use the flux in the normal decay phase to calculate 
Ek.'iso- For the Non sample, no plateau is derived and we 
use any epoch during the normal decay phase to derive 
Ekaso- For GRBs in the Internal sample, there are two 
possibilities. (I) In some cases, a normal decay phase is 
detected after the internal plateau, e.g. GRBs 050724, 
062006, 070724A, 071227, 101219A, and 111121A in Fig¬ 
ure 1. For these bursts, we use the flux at the first data 
point during the normal decay phase to derive £'K,iso- 
(2) For those bursts whose normal decay segment is not 
observed after the rapid decay of the internal plateau at 
later times (the rest of GRBs in Figure 1), we use the 
last data point to place an upper limit to the underlying 
afterglow flux. An upper limit of EK,iso is then derived. 

We adopt two typical values of the circumburst density 
to calculate the afterglow flux, n = 1 cm“^ (a typical 
density of the ISM inside a galaxy) and n = 10“^ cm“^ (a 
typical density in the ISM/IGM with a large offset from 
the galaxy center). For the late epochs we are discussing, 
fast cooling is theoretically disfavored and we stick to 
the slow cooling {vm < Vc) regime. Using the spectral 
and temporal information from the X-ray data, we can 
diagnose the spectral regime of the afterglow based on the 
closure relations (e.g. Zhang & Meszaros 2004; see Gao et 
al. 2013a for a complete review). Most GRBs belong to 
the V > max{vrn,t^c) regime, and we use Equations(lO) 
and (11) of Zhang et al. (2007a) to derive AK,iso- In 
some cases, the spectral regime Vm < v < Vc inferred 
and Equation (13) of Zhang et al. (2007a) is adopted to 
derive EK.iso 

In order to place an upper limit of EK,iso for the Inter¬ 
nal sample GRBs without a detected external shock com¬ 
ponent, one needs to assume the spectral regime and de¬ 
cay slope of the normal decay. To do so, we perform a sta¬ 
tistical analysis of the decay slope and spectral index in 
the normal decay phase using the External and Non sam¬ 
ples (Figure 6). Fitting the distributions with a Gaussian 
distribution, we obtain center values of oq.c = 1.21 ±0.04 
and j3x,c = 0.88±0.05. We adopt these values to perform 
the calculations. Since 2 q!o ~ 3/3jv is roughly satisfied, 
the spectral regime belongs to Vm < v < Vc, and again 
Equation (13) of Zhang et al. (2007a) is again used to 
derive the upper limit of EK.iso- 

In our calculations, the microphysics parameters of the 
shocks are assigned to standard values derived from the 
observations (e.g. Panaitescu & Kumar 2002; Yost et al. 
2003): ee=0.1 and cb = 0.01. The Gompton parameter 
is assigned to a typical value of Y = 1. The calculation 
results are shown in Table 2. 

After obtaining the break time ti, through light curve 
fitting, we derive the bolometric luminosity at the break 


Magnetar central engine in short GRBs 


5 


time tb'. 

Lb = iirDlFb ■ k, (2) 

where Fb is the X-ray flux at tb and k is the /c-correction 
factor. For the Internal sample, we derive the isotropic 
internal plateau energy, Ex,iso^ using the break time and 
break luminosity (Lii & Zhang 2014), i.e. 

Ex,^so ^ Lb ■ (3) 

1 + z 

This energy is also the isotropic emission energy due to 
internal energy dissipation. 

Comparisons of the statistical properties of various de¬ 
rived parameters for the Internal and External samples 
are presented in Figure 7. Figure 7(a) and (b) show 
the distributions of the internal plateau luminosity and 
duration. For the External sample, no internal plateau 
is detected and we place an upper limit on the internal 
plateau luminosity using the observed luminosity of the 
external plateau. The internal plateau luminosity of the 
Internal sample is Lb ergs s“^. The distribution 

of the upper limits of Lb of the External sample peaks 
at a smaller value of Lb lO'^’^-® ergs s“^. This suggests 
that the distribution of internal plateau luminosity Lb 
has an intrinsically very broad distribution (Figure 7(a)). 
The distribution of the duration of the plateaus for the 
Internal sample peaks around 100 s, which is systemat¬ 
ically smaller than the duration of the plateaus in the 
External sample, which peaks around 10^'^ s. In Figure 
7(a) and (b), we also compare the plateau luminosity and 
duration distributions of our sample with those of long 
GRBs (Dainotti et al. 2015) and find that the Internal 
sample is quite different from long GRBs, whereas the 
External sample resembles the distributions of the long 
GRBs well. According to our interpretation, the dura¬ 
tion of the internal plateaus is defined by the collapse 
time of a supra-massive neutron star (Troja et al. 2007; 
Zhang 2014). For the external plateaus, the duration of 
the plateau is related to the minimum of the spin-down 
time and the collapse time of the magnetar. Therefore, 
by definition, the External sample should have a higher 
central value for the plateau duration than the Internal 
sample. The observations are consistent with this hy¬ 
pothesis. 

Figure 7(c) and (d) show the distribution of 7 —ray en¬ 
ergy {E^, iso^ and the internal dissipation energy (^Ex^iso^• 
The E^^iso of the Internal sample is a little bit less than 
that of the External sample, but Ex,iso is much larger 
(for the External sample, only an upper limit of Ex,iso 
can be derived). This means that internal dissipation is 
a dominant energy release channel for the Internal sam¬ 
ple. Figures 7(e) and (f) show the distributions of the 
blastwave kinetic energy (Aic.iso) for different values of 
the number density, n = I cm“^ and n = 10“^ cm“^. 
In both cases, EK,iso of the Internal sample is system¬ 
atically smaller than that for the External sample. The 
results are presented in Tables 2 and 3. 

In Figures 7(g) and (h) (for n = 1,10“^ cm“^, respec¬ 
tively), we compare the inferred total energy of GRBs 
(Atotai = E^ -\- Ex + Ex) with the total rotation energy 
Erot of the millisecond magnetar: 

E,ot = zz 3.5 X 10^2 M2a&RIPo,U^ (4) 


where / is the moment of inertia, R, Pq, and Hq are 
the radius, initial period, and initial angular frequency 
of the neutron star, and M is normalized to the sum 
of the masses of the two NSs ( 2 . 46 M 0 ) in the observed 
NS—NS binaries in our GalaxyQ Hereafter, the conven¬ 
tion Q = lO^Qx is adopted in cgs units for all of the 
parameters except the mass. It is found that the to¬ 
tal energy of the GRBs are below the Ej-ot line if the 
medium density is high (n = 1 cm“^). This energy bud¬ 
get is consistent with the magnetar hypothesis, namely, 
all of the emission energy ultimately comes from the 
spin energy of the magnetar. For a low-density medium 
{n = 10“^ cm“^), however, a fraction of GRBs in the Ex¬ 
ternal sample exceed the total energy budget. The main 
reason for this is that a larger Ex,iso is needed to com¬ 
pensate a small n in order to achieve a same afterglow 
flux. If these GRBs are powered by a magnetar, then the 
data demand a relatively high n. This is consistent with 
the argument that the External sample has a large n and 
so the external shock component is more dominant. 

Figure 8 (a) shows the observed X-ray luminosity at 
t = 10 ^ s (T^^io^s) as a function of the decay slope 02 - 
Figures 8 (b) and 8 (c) show the respective distributions 
of and 02 - The Internal and External samples 

are marked in red and black, respectively. On average, 
the Internal sample has a relatively smaller Lt^igSg than 
that of the External sample (Figure 8 (b)). The fitting 
results from the distributions of various parameters are 
collected in Table 3. 

4. THE MILLISECOND MAGNETAR CENTRAL ENGINE 
MODEL AND IMPLICATIONS 

In this section, we place the short GRB data within 
the framework of the millisecond magnetar central engine 
model and derive the relevant model parameters of the 
magnetar, and discuss the physical implications of these 
results. 

4.1. The millisecond magnetar central engine model 

We first briefly review the millisecond magnetar cen¬ 
tral engine model of short GRBs. After the coalescence 
of the binary NSs, the evolutionary path of the central 
post-merger product depends on the unknown EOS of 
the neutron stars and the mass of the protomagnetar, 
Mp. If Mp is smaller than the non-rotating Tolman- 
Oppenheimer-Volkoff maximum mass MtoVj then the 
magnetar will be stable in equilibrium state (Cook et al. 
1994; Giacomazzo & Perna 2013, Ravi & Lasky 2014). If 
Mp is only slightly larger than MroVi then it may survive 
to form a supra-massive neutron star (e.g. Duez et al. 
2006), which would be supported by centrifugal force for 
an extended period of time, until the star is spun down 
enough so that centrifugal force can no longer support 
the star. At this epoch, the neutron star would collapse 
into a black hole. 

Before the supra-massive neutron star collapses, it 
would spin down due to various torques, the most dom¬ 
inant of which may be the magnetic dipole spin down 

Strictly speaking, M is normalized to the mean of the sum of 
masses of binary NS systems, taking into account the conservation 
of rest mass (Lasky et al. 2014), and ignoring the negligible mass 
lost during the merger process (e.g., Hotokezaka et al. 2013). 
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(Zhang & Meszaros 2001 )E3 The characteristic spin- 
down timescale r and characteristic spin-down luminos¬ 
ity Lq depend on IIq = 2ttIPq and the surface mag¬ 
netic field at the pole Bp, which read (Zhang & Meszaros 
2001 ) 


3c^I 




= 2.05 X 10^ s {hBB-%Pl_^R^^), 


(5) 


Ln — 


2r 


= 1.0 X lO^s erg s-\Bl,^P^X^R%). (6) 


For a millisecond magnetar, the open field line region 
opens a very wide solid angle, and so the magnetar wind 
can be approximated as roughly isotropic. 

Another relevant timescale is the collapse time of a 
supra-massive magnetar, tcoi- For the Internal sample, 
the observed break time R either corresponds to tcoi or 
T, depending on the post-break decay slope a 2 . If 02 — 
2, then the post-break decay is consistent with a dipole 
spin-down model, so that 4 is defined by r and one has 
tcoi > T. On the other hand, if the post-decay slope is 
steeper than 3, i.e. a 2 > S, then one needs to invoke an 
abrupt cessation of the GRB central engine to interpret 
the data (Troja et al. 2007; Rowlinson et al. 2010, 2013; 
Zhang 2014). The break time is then defined by the 
collapse time tcoi, and one has tcoi < t. Overall, one can 
write 


and 


In both cases, the characteristic spin-down luminosity 
is essentially the plateau luminosity, which may be esti¬ 
mated as 


the measured Lq and r can be directly used to derive 
these two parameters. For the Internal sample, both Pq 
and Bp can be derived if a 2 = 2. If a 2 > 3, we can derive 
the upper limit for Pq and The results are presented 
in Table 2 and Figure 9(a)E3 
Figure 9(b) shows the distribution of the collapse times 
for our Internal sample. For GRB 06I20I and GRB 
0707I4B, the decay slope following the plateau is a 2 ^ 2, 
which means that we never see the collapsing feature. A 
lower limit of the collapse time can be set by the last ob¬ 
servational time, so that the stars should be stable long- 
lived magnetars. For the collapsing sample, the center 
value of the tcoi distribution is ~ 100 s, but the half 
width spans about one order of magnitude. 

Figure 10(a) presents an anti-correlation between Lq 
and tcoi, i-e. 

logLo .49 = (-2.79 ± 0.39) logUi .2 - (0.45 ± 0.28) 

( 12 ) 

with r — 0.87 and p < 0.0001. This suggests that a 
longer collapse times tends to have a lower plateau lu¬ 
minosity. This is consistent with the expectation of the 
magnetar central engine model. The total spin energy of 
the millisecond magnetars may be roughly standard. A 
stronger dipole magnetic field tends to power a brighter 
plateau, making the magnetar spin down more quickly, 
and therefore giving rise to a shorter collapse time (see 
also Rowlinson et al. 2014). 

Figure 10(b) presents an anti-correlation between 


= fb/(l + z),a2 = 2, 

> fb/(l + z),a2 > 3. 

(7) 

A^total,iso fmd Gob 

logF;total.iso,52 = (-1.08 ± 0.27) logGol.2 + (0.11 ± 0.18) 

(13) 

f > fb/(l + z),a2 = 2, 

(= fb/(l + z),a2 > 3. 

(8) 

with r = 0.71 and p = 0.0009. This may be understood 
as the follows. A higher plateau luminosity corresponds 
to a shorter spin-down timescale. It is possible that in 


Lq cs Lq 


(9) 


4.2. Magnetar parameters and correlations 

With the above model, one can derive magnetar pa¬ 
rameters and perform their statistics. Two important 
magnetar parameters to define magnetar spin down, i.e. 
the initial spin period Pq and the surface polar cap 
magnetic field Bp, can be solved from the characteristic 
plateau luminosity Lq (Equation (|6])) and the spin-down 
timescale t (Equation JS]); (Zhang & Meszaros 2001), 
i.e. 


Bv,ib = 2.05 G(/45i?g 


-3 


— 1 / 2 - 1 ^ 
^0,49 ^3 ) 


( 10 ) 


^^ 0.-3 = 1 


30 . ^ri/2r-i/2. 


this case, the collapse time is closer to the spin-down 
timescale, and so most energy is already released before 
the magnetar collapses to form a black hole. A lower 
plateau luminosity corresponds to a longer spin-down 
timescale, and it is possible that the collapse time can 
be much shorter than the spin-down timescale, so that 
only a fraction of the total energy is released before the 
collapse. 

Empirically, (Dainotti et al. 2008, 2010, 2013) dis¬ 
covered an anti-correlation between Lf, and tb for long 
GRBs. In Eigure 10(c) we plot our short GRB Internal 
-|- External sample and derive an empirical correlation of 

logLf ,,49 = (-1.41 ± 0.14) log 4.3 - (0.46 ± 0.37), (14) 

with r = 0.88 and p < 0.001. The slope of the correla¬ 
tion is slightly steeper than that of the “Dainotti rela¬ 
tion” (e.g. Dainotti et al. 2008, data see gray dots in 
Fig. 10 (c)). This is probably related to different progen¬ 
itor systems for long and short GRBs, in particular, the 
dominance of Internal plateaus in our sample. Rowlinson 


Since the magnetar wind is likely isotropic for short 
GRBs (in contrast to long GRBs; Lii & Zhang 2014), 

Deviations from the simple dipole spin-down formula may be 
expected (e.g. Metzger et al. 2011; Siegel et al. 2014), but the 
dipole formula may give a reasonable first-order approximation of 
the spin-down law of the nascent magnetar. 


The derived magnetar parameters of most GRBs are slightly 
different from those derived by Rowlinson et al. (2013). One main 
discrepancy is that they used Mp = 1.4M0 to calculate the proto- 
magnetar’s moment of inertia I, wheareas we used Mp = 2.46Mq, 
which is more relevant for post-merger products. The different 
data selection criteria and fitting methods also contribute to the 
discrepancies between the two works. 
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et al. (2014) performed a joint analysis of both long and 
short GRBs taking into account for the intrinsic slope 
of the luminosity-time correlation (Dainotti et al. 2013). 
We focus on short GRBs only but studied the Internal 
and External sub-samples separately. 


4.3. Constraining the neutron star EOS 

The inferred collapsing time can be used to constrain 
the neutron star EOS (Lasky et al. 2014; Ravi & Lasky 
2014). The basic formalism is as follows. 

The standard dipole spin-down formula gives (Shapiro 
& Teukolsky 1983) 


^W = Po(l + ^^t)V^ 


3c3 
t , 


= Po(l + -)^/" 


(15) 


For a given EOS, a maximum NS mass for a non-rotating 
NS, i.e. MtoV) can be derived. When an NS is supra- 
massive but rapidly rotating, a higher mass can be sus¬ 
tained. The maximum gravitational mass (Mmax) de¬ 
pends on spin period, which can be approximated as 
(Lyford et al. 2003) 


Internal sample with redshift measurements. Five NS 
equations of state, i.e. SLy (black, Douchin & Haensel. 
2001), APR (red, Akmal et al. 1998), GMl (green. Glen- 
denning & Moszkowski. 1991), AB-N, and AB-L (blue 
and cyan, Arnett & Bowers. 1997) are shown in different 
vertical color bands. The gray shaded region is the pro- 
tomagnetar mass distribution, Mp, discussed above. The 
horizontal dashed line is the observed collapse time for 
each short GRB. Our results show that the GMl model 
gives an Mp band fall in the 2cr region of the protomagne- 
tar mass distribution, so that the correct EOS should be 
close to this model. The maximum mass for non-rotating 
NS in this model is Mrov = 2.37Mq. 

Lasky et al. (2014) applied the observational collapse 
time of short GRBs to constrain the NS EOS (see also 
a rough treatment by Fan et al. 2013a). Our results 
are consistent with those of Lasky et al. (2014), but 
using a larger sample. Another improvement is that we 
introduce a range of t rather than one single r to derive 
the range of plausible Mp, since the observed collapse 
time only gives the lower limit of r. This gives a range 
of the allowed Mp (rather than a hne-tuned value for the 
single r scenario) for each GRB for a given observed tb- 

5. CONCLUSIONS AND DISCUSSION 


M„,, = MTov(l + a-P'') (16) 


where a and /3 depend on the EOS. The numerical values 
of a and /3 for various EOSs have been worked out by 
Lasky et al. (2014), and are presented in Table 4 along 
with MtoV; hS, and I. 

As the neutron star spins down, the maximum mass 
Afmax gradually decreases. When Mmax becomes equal 
to the total gravitational mass of the protomagnetar, Mp, 
the centrifugal force can no longer sustain the star, and 
so the NS will collapse into a black hole. Using equa¬ 
tion Equation (US]) and Equation (ITOl) . one can derive the 
collapse time: 


^col- 


Sc^J r,Mp — Mtov 


47r2p2^6 


[(- 


AMtov 


)2//i _ 


_ u Mp-MTOY ^2/0 _ p2i 
Pq cIMtov ° 


(17) 


As noted, one can infer Bp, Pq, and Uoi from the obser¬ 
vations. Moreover, as the Galactic binary NS population 
has a tight mass distribution (e.g., Valentim et al. 2011; 
Kiziltan et al. 2013), one can infer the expected distri¬ 
bution of protomagnetar masses, which is found to be 
Mp = 2A6^1 \^Mq (for details see Lasky et al. 2014). 
The only remaining variables in Equation (16) are re¬ 
lated to the EOS, implying that the observations can be 
used to derive constraints on the EOS of nuclear matter. 
For most GRBs in our Internal sample, only the lower 
limit of T is derived from 4 (Equation ([T])). One can 
also infer the maximum r by limiting Pq to the break-up 
limit. Considering the uncertainties related to gravita¬ 
tional wave radiation, we adopt a rough limit of 1 mil¬ 
lisecond. By doing so, one can then derive a range of r, 
and hence a range of Mp based on the data and a given 
EOS. 

Figure 11 presents the collapse time (Uoi) as a function 
of protomagnetar mass {Mp) for each short GRB in the 


In this paper, by systematically analyzing the BAT- 
XRT light curves of short GRBs detected by Swift be¬ 
fore 2014 August, we systematically examine the mil¬ 
lisecond magnetar central engine model of short GRBs. 
About 40 GRBs have bright X-ray afterglows detected 
with Swift/XBT, 8 Of which have EE detected with 
Swift/BAT. Based on the existence of plateaus, their ob¬ 
servation properties, and how likely a GRB is powered 
by a millisecond magnetar central engine, we character¬ 
ized short GRBs into three samples: Internal (plateau), 
External (plateau), and Non (plateau). We compared 
the statistical properties of our samples and derived or 
placed limits on the magnetar parameters Pq and Bp 
from the data. Using the collapse time Uoi of the proto¬ 
magnetar inferred from the plateau break time tb in the 
Internal sample, we went on to constrain the NS EOS. 
The following interesting results are obtained. 

• At least for the Internal sample, the data seem to 
be consistent with the expectations of the mag¬ 
netar central engine model. Assuming isotropic 
emission, the derived magnetar parameters Bp and 
Pq fall into a reasonable range. The total energy 
(sum of E-y, Ex, and Ek) is within the budget pro¬ 
vided by the spin energy of the millisecond mag¬ 
netar {Erot ~ 3-5 X 10®^ erg). The Lq — tcoi anti¬ 
correlation is generally consistent with the hypoth¬ 
esis that the total spin energy of the magnetar may 
be standard, and a higher dipolar magnetic field 
powers a brighter but shorter plateau. 

• The so-called EE following some short GRBs is es¬ 
sentially the brightest internal plateau commonly 
observed in short GRBs. A more sensitive and 
softer detector would detect more EE from short 
GRBs. 

• The External sample may also be consistent with 
having a magnetar central engine, even though the 
evidence is not as strong. If both the Internal 
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and External samples are powered by a millisec¬ 
ond magnetar central engine, then the difference 
between the two samples may be related to the cir- 
cumburst medium density. The physical and host- 
normalized offsets of the afterglow locations for the 
Internal sample is somewhat larger than those of 
the External sample, even though the separation 
between the two samples is not clear cut. In any 
case, it is consistent with this expectation. The 
total energy budget of the GRB is within the mag¬ 
netar energy budget for the External sample only 
if the ambient density is relatively large, and hence 
powers a strong external shock emission compo¬ 
nent. There is no significant difference between 
those two groups for the star formation rate, metal- 
licity, and age of the host galaxy. 

• Using the collapse time of supra-massive protomag- 
netar to form a black hole and the distribution of 
the total mass of NS—NS binaries in the Galaxy, 
one can constrain the NS EOS. The data point to¬ 
ward an EOS model close to GMI, which has a non¬ 
spinning maximum NS mass Mtov ~ 2.37Mq. 

The short GRB data are consistent with the hypoth¬ 
esis that the post-merger product of an NS—NS merger 
is a supra-massive neutron star. The existence of such 
a long-lived post-merger product opens some interest¬ 
ing prospects in the multi-messenger era. In particular, 
the dipole spin-down power of the supra-massive NS can 
power bright electromagnetic radiation even if the short 
GRB jet does not beam toward earth, and so some in¬ 
teresting observational signatures are expected to be as¬ 
sociated with gravitational wave signals in the Advanced 
LIGO/Virgo era (Fan et al. 2013b; Gao et al. 2013; Yu 
et al. 2013; Zhang 2013; Metzger & Piro 2014). An¬ 
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other interesting possibility is that a fast radio burst 
(e.g. Lorimer et al. 2007; Thornton et al. 2013) may 
be released when the supra-massive magnetar collapses 
into a black hole (Falcke & Rezzolla 2014; Zhang 2014). 
The discovery of an FRB following a GRB at the end 
of the internal plateau (see Bannister et al. 2012) would 
nail down the origin of FRBs, although such observations 
require fast telescope response times given the expected 
distribution of collapse times following SGRBs (see figure 
9(b) and Ravi & Lasky 2014). The GRB-FRB associa¬ 
tions, if proven true, would be invaluable for cosmology 
studies (Deng & Zhang 2014; Gao et al. 2014; Zheng et 
al. 2014 Zhou et al. 2014). 

Recently, Rezzolla & Kumar (2014) and Giolfi & Siegel 
(2014) proposed a different model to interpret the short 
GRB phenomenology. In their model, the post-merger 
product is also a supra-massive NS, but the collapse time 
is allocated as the epoch of the short GRB itself, rather 
than the end of the Internal plateau. Our conclusions 
drawn in this paper do not apply to that model. A crucial 
observational test to differentiate between our model and 
theirs is whether or not there exists strong X-ray emission 
before the short GRB itself. This may be tested in the 
future with a sensitive wide-field XRT. 
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TABLE 1 

Observed properties of short GRBs in our samples. 


GRB 

name 


Tgo/EEb 

(s) 

r ^ 

i ry 


Host Offset® 
(kpc) 

Host Offset® 

(Te) 

4 ' 

(s) 

0 . 1 ^ 

0 . 2 ^ 

X^/dof 

Ref 

Internal 

050724 

0.2576 

3/154 

1.89±0.22 

0.58±0.19 

2.76±0.024 

— 

139±9 

0.20±0.1 

4.16±0.05 

980/835 

(1, 

051210 

(0.58) 

1.27/40 

1.06±0.28 

1.1±0.18 

24.9±24.6 

4.65±4.6 

67±4 

0.15±0.04 

2.96±0.09 

118/132 

(1, 2) 

051227 

(0.58) 

3.5/110 

1.45±0.24 

1.1±0.4 

— 

— 

89±5 

0.10±0.05 

3.19±0.13 

681/522 

(3, 4, 5 

060801 

1.13 

0.49/N 

1.27±0.16 

0.43±0.12 

— 

— 

212±11 

O.lOiO.ll 

4.35±0.26 

81/75 

(1) 

061006 

0.4377 

0.5/120 

1.72±0.17 

0.76±0.28 

1.3±0.24 

0.35±0.07 

99±7 

0.17±0.03 

9.45±1.14 

111/138 

(1> 2) 

061201 

0.111 

0.76/N 

0.81±0.15 

1.2±0.22 

32.47±0.06 

14.91±0.03 

2223±43 

0.54±0.06 

1.84±0.08 

20/24 

(1> 6) 

070714B 

0.9224 

3/100 

1.36±0.19 

1.01±0.16 

12.21±0.53 

5.55±0.24 

82±2 

0.10±0.07 

1.91±0.03 

672/581 

(1> 6) 

070724A 

0.46 

0.4/N 

1.81±0.33 

0.5±0.3 

5.46±0.14 

1.5±0.04 

77±6 

O.OliO.l 

6.45±0.46 

256/222 

(1, 6) 

071227 

0.381 

1.8/100 

0.99±0.22 

0.8±0.3 

15.5±0.24 

3.28±0.05 

69±8 

0.27±0.08 

2.92±0.06 

244/212 

(1, 6) 

080702A 

(0.58) 

0.5/N 

1.34±0.42 

1.03±0.35 

— 

— 

586±14 

0.51±0.22 

3.56±0.31 

3/5 

(7) 

080905A 

0.122 

1/N 

0.85±0.24 

0.45±0.14 

17.96±0.19 

10.36±0.1 

13±3 

0.19±0.09 

2.37±0.07 

43/52 

(6, 7) 

080919 

(0.58) 

0.6/N 

1.11±0.26 

1.09±0.36 

— 

— 

340±26 

0.40±0.14 

5.20±0.55 

7/5 

(7) 

081024A 

(0.58) 

1.8/N 

1.23±0.21 

0.85±0.3 

— 

— 

102±5 

0.27±0.02 

5.89±0.3 

50/42 

(7) 

090510 

0.903 

0.3/N 

0.98±0.21 

0.75±0.12 

10.37±2.89 

1.99±0.39 

1494±87 

0.69±0.04 

2.33±0.11 

112/132 

(6, 7) 

090515 

(0.58) 

0.036/N 

1.61±0.22 

0.75±0.12 

75.03±0.15 

15.53±0.03 

178±3 

0.10±0.08 

12.62±0.5 

42/38 

(6, 7) 

100117A 

0.92 

0.3/N 

0.88±0.22 

1.1±0.26 

1.32±0.33 

0.57±0.13 

252±9 

0.55±0.03 

4.59±0.13 

84/92 

(6, 7, 8 

100625A 

0.425 

0.33/N 

0.91±0.11 

1.3±0.3 

— 

— 

200±41 

0.26±0.44 

2.47±0.18 

3/6 

(7, 9) 

100702A 

(0.58) 

0.16/N 

1.54±0.15 

0.88±0.11 

— 

— 

201±6 

0.62±0.13 

5.28±0.23 

82/69 

(7) 

101219A 

0.718 

0.6/N 

0.63±0.09 

0.53±0.26 

— 

— 

197±10 

0.13±0.19 

20.52±8.01 

3/5 

(7) 

111121A 

(0.58) 

0.47/119 

1 .66±0.12 

0.75±0.2 

— 

— 

56±9 

0.10±0.13 

2.26±0.04 

274/289 

(7) 

120305A 

(0.58) 

0.1/N 

1.05±0.09 

1.4±0.3 

— 

— 

188±14 

0.73±0.14 

6.49±0.63 

14/18 

(7) 

120521A 

(0.58) 

0.45/N 

0.98±0.22 

0.73±0.19 

— 

— 

270±55 

0.30±0.27 

10.74±4.76 

3/7 

(7) 

External 

051221A 

0.55 

1.4/N 

1.39±0.06 

1.07±0.13 

1.92±0.18 

0.88±0.08 

25166±870 

0.12±0.13 

1.43±0.04 

52/63 

(1, 2, 7 

060313 

(0.58) 

0.71/N 

0.71±0.07 

1.06±0.15 

2.28±0.5 

1.23±0.23 

2294±65 

0.3±0.15 

1.52±0.04 

54/45 

(1, 2, 7 

060614 

0.1254 

5/106 

2.02±0.04 

1.18±0.09 

— 

— 

49840±3620 

0.18±0.06 

1.9±0.07 

70/54 

(1,10) 

070714A 

(0.58) 

2/N 

2 .6±0.2 

1.1±0.3 

— 

— 

892±34 

0.11±0.09 

0.95±0.06 

15/18 

(7) 

070809 

0.219 

1.3/N 

1.69±0.22 

0.37±0.21 

33.22±2.71 

9.25±0.75 

8272±221 

0.18±0.06 

1.31±0.17 

17/22 

(6, 7) 

080426 

(0.58) 

1.7/N 

1.98±0.13 

0.92±0.24 

— 

— 

566±97 

0.11±0.16 

1.29±0.05 

28/21 

(7) 

090426 

2.6 

1.2/N 

1.93±0.22 

1.04±0.15 

0.45±0.25 

0.29±0.14 

208±53 

0.12±0.07 

1.04±0.04 

15/11 

(6, 7) 

100724A 

1.288 

1.4/N 

1.92±0.21 

0.94±0.23 

— 

— 

5377±331 

0.72±0.08 

1.61±0.12 

16/19 

(11, 12' 

130603B 

0.356 

0.18/N 

1.83±0.12 

1.18±0.18 

5.21±0.17 

1.05±0.04 

3108±356 

0.4±0.02 

1.69±0.04 

126/109 

(6, 13, 

130912A 

(0.58) 

0.28/N 

1 .21±0.2 

0.56±0.11 



231±54 

0.04±0.39 

1.34±0.04 

28/21 

(15) 


References. — (l)Zhang et al.(2009),(2)Fong et al.(2010),(3)Hullinger et al.(2005),(4)Butler et al. (2007),(5)Gompertz et al. (2014),(6)Pong Sc 

Berger(2013), (7)Rowlinson et al. (2013),(8)Pong et al. (2011), (9)Pong et al. (2013),(10)Lu Sc Zhang. (2014), (1 l)Tlioene et al. (2010), (12)Markwardt et 
al.(2010),(13)Cucchiara et al.(2013),(14)Barthelmy et al.(2013),(15)Krimm et al.(2013). 


^ The measured redshift are from the published papers and GNCs. When the redshift is not known, 0.58 is used. 

The duration of the GRB without and with extended emission (if EE exists). “N” denotes no EE. 

^ The photon index in the BAT band (15-150keV) fitted using a power law. 

The spectral index of the absorbed power-law model for the normal segments. 

^ Physical and host-normalized offsets for the short GRBs with Hubble Space Telescope (HST) observations, 
f 

The break time of the light curves from our fitting, and Q 2 the decay slopes before and after the break time. 
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TABLE 2 

The derived properties oe the short GRBs in our samples. 


GRB 

Name 

iso, 51 
(erg) 

r b 

-^6,49 
(erg s-l) 

(S) 

Bp,15^ 

(G) 

4),-3" 
(s-1) 

i47(103s)d 
(erg s-l) 

^K,iso,51^ 
(n = l,erg) 

^K,iso,51^ 

(n = 10“^, erg) 

^X,iso,51^ 
(erg) 

Internal 

050724 

u.uy_Q Q2 

l.litO.16 

o.llt 

17.15t 

4.044. 

0.05±0.006 

0.97±0.13 

2.37±0.26 

1.25±0.14 

051210 

r» Qr)-|-0.036 
U.^^_0.036 

2.23±0.26 

0.04t 

32.69t 

4.684. 

0.324. 

0.34t 

1.894. 

0.94±0.10 

051227 


1.89±0.02 

0.05t 

28.68t 

4.574. 

0.66±0.086 

2.69±0.35 

5.65±0.26 

0.98±0.11 

060801 

1.70t°'i 

0.73±0.07 

0.14t 

17.8U 

4.57i 

0.46i 

3.84t 

2.034. 

0.98±0.16 

061006 

2.20t\l 

3.37±0.32 

0.04t 

18.3U 

3.164. 

0.17±0.022 

6.37±0.83 

6.37±0.83 

2.06±0.23 

061201 

n 1 q+O-02 

(l±0.11)E-3 

2±0.043 

31.17±2.36 

30.80±1.97 

0.15±0.019 

0.74±0.10 

1.84±0.21 

0.02±0.01 

070714B 

usotti 

6.22±0.09 

0.04±0.002 

19.12±1.08 

2.77±0.09 

1.40±0.182 

4.40±0.57 

9.47±0.41 

2.67±0.29 

070724A 

0.03l°;°l 

13.1±7.2 

0.05t 

10.89t 

1.73i 

0.05±0.007 

7.99±1.04 

7.99±1.04 

6.81±4.56 

071227 

2.20l°;| 

0.77±0.01 

0.05t 

44.93t 

7.164. 

0.05±0.007 

0.91±0.12 

2.07±0.23 

0.40±0.04 

080702A 

n 1 O+0.208 

U.10_o.o556 

(7±0.25)E-3 

0.37t 

64.44 

27.40t 

0.02±0.002 

0.26t 

0.84i 

0.03±0.01 

080905A 

7+11 

‘ -4 

2.76±0.9 

o.oit 

102.834. 

7.874. 

0.04 

0.37t 

0.724. 

0.33±0.18 

080919 

^•^^-0.278 

0.05±0.01 

0.22t 

41.98t 

13.63t 

0.11±0.014 

0.20t 

1.04 

0.11±0.04 

081024A 

u.oe_o 278 

0.78±0.14 

0.06t 

36.34t 

6.424. 

0.44 

0.50t 

0.954. 

0.50±0.12 

090510 

0+5 

’^-2 

0.18±0.03 

0.75t 

6.364. 

3.854. 

7.90±1.027 

3.71±0.48 

7.79±0.85 

1.38±0.37 

090515 

0.08l°;J®2 

1.24±0.05 

O.llt 

16.29t 

3.824. 

0.284. 

0.83t 

0.934. 

1.40±0.10 

100117A 

2.50t°;| 

0.45±0.04 

0.16t 

19.06t 

5.32i 

0.02i 

0.12t 

0.92i 

0.72±0.09 

100625A 

U.DA_o 031 

0.042±0.03 

0.13t 

79.67t 

19.75t 

0.02±0.003 

0.07t 

0.14 

0.05±0.05 

100702A 

0.47l°;°« 

0.97±0.14 

0.13t 

16.36t 

4.07i 

1.20i 

1.82t 

4.041 

1.24±0.22 

101219A 

4.80l°;| 

0.56±0.05 

0.12t 

23.84t 

5.654. 

0.234. 

4.14±0.54 

10.03±1.11 

0.64±0.10 

111121A 

2.80t°'i« 

14±0.8 

0.04t 

15.22t 

2.021 

1.57±0.204 

9.80±1.27 

22.64±1.49 

5.04±0.55 

120305A 

n qq+0.0112 

u.^y_o.oii2 

0.48±0.09 

0.13t 

23.30t 

5.8O4. 

0.11±0.014 

0.45±0.06 

0.89±0.10 

0.61±0.17 

120521A 

o qq+ 0.0115 
U.^o_0 0356 

0.07±0.003 

0.17t 

44.68t 

12.90t 

1.04 

2.46t 

4.42i 

0.12±0.03 

External 

051221A 

2.80tf i 

(1.78±0.09)E-5 

— 

— 

— 

0.63±0.08 

16.29±2.12 

35.56±3.91 

0.31±0.032 

060313 

12.90l?|f 

(2.74±0.21)E-2 

— 

— 

— 

3.00±0.39 

8.11±1.05 

17.21±1.89 

0.45±0.054 

060614 

2.40l°'l 

(2.55±0.12)E-4 

— 

— 

— 

0.04±0.01 

7.06±0.92 

14.56±1.61 

0.11±0.012 

070714A 


(1.3±0.15)E-2 

— 

— 

— 

0.70±0.09 

13.94±1.81 

13.94±1.81 

0.07±0.013 

070809 

0 01 

U.Ui_o.oi 

(3.2±0.31)E-5 

— 

— 

— 

0.05±0.01 

2.25±0.29 

5.61±0.62 

0.02±0.003 

080426 

n Q9+1-25 
U.O^_0.0556 

(3.53±1.01)E-2 

— 

— 

— 

0.83±0.11 

4.71±0.61 

10.35±1.14 

0.13±0.064 

090426 

4-201®,4 

2.46±0.48 

— 

— 

— 

12.50±1.63 

52.84±6.87 

128.09±14.09 

1.43±0.690 

100724A 

0.7t°i 

(2.85±0.32)E-2 

— 

— 

— 

5.20±0.68 

13.83±1.80 

30.42±3.34 

0.67±0.180 

130603B 

2.201°'2 

(1.2±0.05)E-2 

— 

— 

— 

1.60±0.21 

6.12±0.80 

13.46±1.26 

0.27±0.055 

130912A 

p 

CO 

1 + 
o o 
b b 
00 00 

0.21±0.09 

— 

— 

— 

1.50±0.20 

20.15±2.62 

49.53±5.45 

0.30±0.237 


^ ^gQ is calculated using fluence and redshift extrapolated into l-10,000keV (rest frame) with a spectral model and a fc-correction, in units of 10^^ erg. 

^ Isotropic luminosity at the break time (in units of 10^® erg s~^) and the spin-down time (in units of lO^s). 

^ The dipolar magnetic field strength at the polar cap in units of lO^^G, and the initial spin period of the magnetar in units of milliseconds, with an assumption of an 


isotropic wind. 

The luminosity of the afterglow at t = 1000 s. The arrow sign indicates the upper limit. 

^ The isotropic kinetic energy measured from the afterglow flux during the normal decay phase, in units of 10^^ erg. 
^ The isotropic internal dissipation energy in the X-ray band (also internal plateau), in units of 10^^ erg. 


TABLE 3 

The center values and standard deviations of the Gaussian fits of various distributions. 


Name 

Internal 

External 

log(L6) erg s-l 

(49.06 ± 0.15) erg s"! 

(47.55 ±0.16) erg s-l 

logith) s 

(2.01 ±0.06) s 

(3.41 ± 0.04) s 

log{Ej, iso) erg 

log{E^, iso, n = 1 cm-3) erg 

(50.78 ±0.16) erg 

(51.25 ±0.08) erg 

(50.86 ± 0.11) erg 

(51.35 ±0.04) erg 

log{Ex, iBo,n = lO”® cm-3) erg 

(51.74 ±0.18) erg 

(52.32 ± 0.06) erg 

log{Etotal,iso,n = 1 cm-3) erg 

(51.36 ± 0.06) erg 

(51.82 ± 0.04) erg 

log{Etota,l,iso,n = 10“® cm-^) erg 

(51.61 ± 0.07) erg 

(52.39 ± 0.03) erg 

log{too\) s 

(1.96 ±0.02) s 

— 

*off(it^l03 J erg s-l 

(46.09 ± 0.07) erg s"! 

(47.08 ±0.09) erg s"! 












12 


Lii et al 


TABLE 4 

The parameters of various NS EOS models 



SLy 

APR 

GMl 

AB-N 

AB-L 

Mrovi^e) 

2.05 

2.20 

2.37 

2.67 

2.71 

R (km) 

9.99 

10.0 

12.05 

12.9 

13.7 

7(10^® g cm^) 

1.91 

2.13 

3.33 

4.30 

4.70 

1 

O 

1 

o 

1.60 

0.303 

1.58 

0.112 

2.92 

/3 

-2.75 

-2.95 

-2.84 

-3.22 

-2.82 


References. — The neutron star EOS parameters are derived in Lasky et al. (2014) and Ravi & Lasky (2014). 
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Fig. 1. — The BAT-XRT rest-frame light curves of the GRBs in our Internal sample. Blue triangles are BAT data extrapolated to the 
XRT band, and black points (with error bars) are the XRT data. The red solid curves are the best fits with a smooth power-law model to 
the data. The green dotted lines are the best fits with power-law model after the steeper decay. 
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Fig.l— continued. 
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Fig. 2.— Similar to Figure 1, but for the External sample, 
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Fig. 3. — Ensemble of X-ray light curves (0.3-10 keV) of the GRBs in our Internal sample with EE, Internal sample without EE, External 
sample, and Non sample. 
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Fig. 4.— (a) Duration distributions of the extended emission for EE sample, and the internal plateau emission for the no-EE sample. 
Inset: the cumulative duration distributions for the EE and no-EE sub-samples, (b) The redshift distribution of all short GRBs with 2 : 
measurements. The red solid line is the best Gaussian fit with a center value Zc = 0.58. The green and blue histograms are the redshift 
distributions for the EE and no-EE sub-samples, respectively, (c) The plateau flux distributions of both EE (red, solid line -|- bar) and 
no-EE (black, solid line) GRBs in our Internal sample. The dotted lines are the best Gaussian fits to the distributions, (d) A joint fit to 
the flux distribution of all the GRBs in the Internal plateau (both EE and no-EE included), (e) The plateau luminosity distributions of 
both EE (red, solid line -|- bar) and no-EE (black, solid line) GRBs in our Internal sample. The dotted lines are the best Gaussian fits to 
the distributions, (f) A joint fit to the luminosity distribution of all the GRBs in the Internal plateau (both EE and no-EE included). The 
insets in (e) and (f) are for the GRBs with measured redshifts only. 
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Fig. 5.— Distributions of the physical offsets and host-normalized offsets of the Internal and External samples. The solid and dash lines 
are the best Gaussian fitting for the Internal and External samples, respectively. 
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Fig. 6.— Distributions of the decay slope ao and spectral index /3x in the normal decay phase in our External and Non samples. The 
solid lines are the best Gaussian fits to the distributions. 
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Fig. 7.— Comparisons of various properties between the Internal (red, open histogram) sample and the External (black, gray histogram) 
samples. The best-fit Gaussian profiles are over-plotted with the respective colors. The eight panels denote histograms of L^, tb, E^^iso^ 
Ex,iso^ ^K,iso nnd Etotal,iso 5 respetively, with the last two parameters plotted twice for two different medium densities, n = 1,10“^ cm~^. 
The vertical dotted line in panels (g) and (h) denotes the total rotation energy budget of a millisecond magnetar. If no redshift is measured, 
then 2 = 0.58 is adopted in the calculations. 
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Fig. 8. — One-dimensional (panels (b) and (c)) and two-dimensional (panel (a)) L(t = 1000 s) — a 2 distributions of the GRBs in our 
samples. The red diamonds and black dots denote the Internal and External samples, respectively, and the arrows indicate the upper limits. 
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Fig. 9. — (a) Inferred magnetar parameters, initial spin period Pq vs. surface polar cap magnetic field strength Bp derived for our Internal 
sample. The red diamonds indicate GRB 061201 and GRB 070714B, which have r measured from tf,. All the other GRBs only have the 
lower limit of r. The arrows denote upper limits. The vertical solid line is the break-up spin period limit for a neutron star (Lattimer Sc 
Prakash 2004). (b) The distribution of the collapse time for our Internal sample. The dotted line is the best Gaussian profile fit. 
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Fig. 10.— (a) Lo — tcoi anti-correlation for our Internal samples. The red diamonds are GRB 061201 and GRB 070714B, and the arrows 
denote the lower limits of the collapse time, (b) The -Ftotal iso ~ ^col anti-correlation for our Internal sample using n = 1 cm“^ to calculate 
^K,iso- The blue diamonds indicate the upper limits to calculate -FK,iso? the red solid line is the best-fitting line, (c) The empirical 
Lfe —tfj correlation derived from the short GRBs in our sample (red for External and blue for Internal samples) compared with the Dainotti 
relation for long GRBs (gray). The solid line is the best power-law fit to the SGRBs sample, and the two dotted lines denote the 2cr region 
of the fit. 
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Fig. 11. — Collapse time as a function of the protomagnetar mass. The shaded region is the protomagnetar mass distribution derived 
from the total mass distribution of the Galactic NS—NS binary systems. The predicted results for 5 equations of state are shown in each 
panel: SLy (black), APR (red), GMl (green), AB-N (blue), and AB-L (cyan). The horizontal dotted line is the observed collapse time for 
each GRB. 




















































